(*) 

Measurement of 5 ^ Dl Ktt branching fractions 



J. Wiechczynski,^'' T. Lesiak,^'''^ H. Aihara,^° K. Arinstein/ V. Aulchenko/ 
A. M. Bakich,^^ V. Balagura/^ E. Barberio,^^ A. Bay/*^ K. Belous/^ V. Bhardwaj,^^ 
A. Bondar/ A. Bozek,^^ M. Bracko,^^'!^ j Brodzicka,* T. E. Browder/ Y. Chao,^^ 

A. Chen,24 B. G. Cheon," R. Chistov/^ I._S. Cho,^^ Y. Choi,^^ J. Dalseno,^ M. Dash,^^ 
A.. Drutskoy,'^ S. Eidelman/ D. Epifanov/ N. Gabyshev/ A. Garmash,^^ P. Goldenzweig; 

H. Ha/6 B.-Y. Han/6 Y. Hoshi,^^ W.-S. Hou,^^ H. J. Hyun/^ T. lijima,^^ K. Inami,22 
A. Ishikawa,33 M. Iwasaki,^° D. H. Kah/^ J. H. Kang,^^ P. Kapusta,^^ H. Kawai,^ 
T. Kawasaki,29 H. Kichimi,^ H. O. Kim/^ J. H. Kim,^^ Y. I. Kim/^ Y. J. Kim,^ 

B. R. P. Krizan/9'14 P. Krokovny,*^ A. Kuzmin/ Y.-J. Kwon/^ S.-H. Kyeong,^^ 
S. E. Lee,35 C. Liu,^^ Y. Liu,^^ D. Liventsev/^ R. Louvot/^ A. Matyja,^^ S. McOnie,^^ 

T. Medvedeva/3 H. Miyata,^^ Y. Miyazaki,^^ T. Mori,^^ Y. Nagasaka,^ E. Nakano,^^ 
M. Nakao,8 H. Nakazawa,^^ K. Nishimura/ O. Nitoh,^^ Ogawa,^^ T. Ohshima,^^ 
S. Okuno/5 H. Ozaki,8 P. Pakhlov/^ G. Pakhlova/^ H. Palka,^^ C. W. Park,^^ 

H. K. Park/7 K. S. Park,^^ R. Pestotnik/^ L. E. Piilonen,^^ Poluektov/ 
H. Sahoo/ Y. Sakai,8 O. Schneider/^ C. Schwanda/^ A. Sekiya,23 K. Senyo,22 
M. E. Sevior,2i M. Shapkin/^ v. Shebalin/ J.-G. Shiu,^^ B. Shwartz/ J. B. Singh,32 
A. Sokolov/2 S. Stanic,30 M. Staric/^ J. Stypula,^^ T. Sumiyoshi,^^ G. N. Taylor,^! 
Y. Teramoto,'^^ I. Tikhomirov/^ S. Uehara,^ K. Ueno,^^ T. Uglov/^ Y. Unno,^ S. Uno,* 

Y. Usov/ G. Varner/ K. E. Varvell,^'' K. Vervink/^ A. Vinokurova/ C. H. Wang,^^ 
M.-Z. Wang,26 R Wang/° Y. Watanabe/^ E. Won/^ B. D. Yabsley,^^ Y. Yamashita,^^ 
V. Zhilich/ V. Zhulanov/ T. Zivko/^ A. Zupanc/'* N. Zwahlen/^ and O. Zyukova^ 

(The Belle Collaboration) 

^Budker Institute of Nuclear Physics, Novosibirsk 

^Chiba University, Chiba 
^University of Cincinnati, Cincinnati, Ohio 45221 
^ T. Kosciuszko Cracow University of Technology, Krakow 
^The Graduate University for Advanced Studies, Hayama 
^Hanyang University, Seoul 
'^University of Hawaii, Honolulu, Hawaii 96822 
^High Energy Accelerator Research Organization (KEK), Tsukuba 
'^Hiroshima Institute of Technology, Hiroshima 
^^Institute of High Energy Physics, 
Chinese Academy of Sciences, Beijing 
^^Institute of High Energy Physics, Vienna 
^^Institute of High Energy Physics, Protvino 
^^Institute for Theoretical and Experimental Physics, Moscow 
^^J. Stefan Institute, Ljubljana 
^^Kanagawa University, Yokohama 
i^Xorea University, Seoul 
^'^Kyungpook National University, Taegu 
^^Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne 



Typeset by REVTgX 



1 



^^Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana 

University of Maribor, Maribor 
University of Melbourne, School of Physics, Victoria 3010 
^^Nagoya University, Nagoya 
^^Nara Women's University, Nara 
^"^National Central University, Chung-li 
^^National United University, Miao Li 
Department of Physics, National Taiwan University, Taipei 
^'^H. Niewodniczanski Institute of Nuclear Physics, Krakow 
^^Nippon Dental University, Niigata 
^^Niigata University, Niigata 
University of Nova Gorica, Nova Gorica 
Osaka City University, Osaka 
^^Panjab University, Chandigarh 
^^Saga University, Saga 
^'^ University of Science and Technology of China, Hefei 
^^Seoul National University, Seoul 
Sungkyunkwan University, Suwon 
^'''University of Sydney, Sydney, New South Wales 
^^Toho University, Funabashi 
^^Tohoku Cakuin University, Tagajo 
'^'^ Department of Physics, University of Tokyo, Tokyo 

Tokyo Metropolitan University, Tokyo 
Tokyo University of Agriculture and Technology, Tokyo 
^^IPNAS, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24.061 

^'^ Yonsei University, Seoul 
Princeton University, Princeton, New Jersey 08544 

Abstract 

We report a measurement of the exclusive meson decay to the oi*^ K+Tr"*" final state ris- 
ing 657 X 10^ BB pairs collected at the T(4S') resonance with the Belle detector at the KEKB 
asymmetric-energy e"^e~ collider. We use D*^ -^77 ^7 ~^ j K*{892)^K~ 

and K^K~ decay modes for D^*^ reconstruction and measure the following branching fractions: 
B{B+ D~K+Tr+) = (IMtofsi^tat) tofoisyst) ± 0.17(B)) x 10"^ and B{B+ D*-K+Tr+) = 
{lA7'^Q']^{stat) iQjg(syst) lb 0.13(B)) x 10~^. The uncertainties are due to statistics, experimental 
systematic errors and uncertainties of intermediate branching fractions, respectively. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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The dominant process in the decays Ds* K^n^ [1] is mediated by the b ^ c quark 

transition and includes the production of an additional ss pair, as shown by the Feynman 
diagram in Fig. 1(a). This process produces at least three final-state particles and can thus be 
distinguished from much more dominant decays, which include direct Ds production from the 
W boson vertex. An example of a process that does not involve ss quark popping is shown in 
Fig. 1(b); this is the dominant Feynman diagram describing two-body D^*^^ D decays 

with D° ^ K+TT~. Although both B+ Di*^^K+TT+ and B+ Di*^^D^(D° K+tt~) 
decays give a similar three-body final state, the different decay mechanisms lead to opposite 
charges for the Dg and vr mesons. In addition, due to the similarities of the final states, 
the latter decay, B+ Di*^'^D^, can be used to check the experimental procedure for the 
exclusive measurements of the former one, B^ Di*^~ K'^tv^ . These three body decay 
modes were recently observed by BaBar [2] and need further confirmation. 

Studies of B~^ — > Di*^~ K^ir^ decays are additionally motivated by interest in the inter- 
mediate resonances that can be formed from the three final-state particles. These resonances 
would be visible in the Dalitz plots for different two-body subsystems [3]. 

In this paper we report measurements of the branching fractions for B^ oi*^ K^tx^ 
decays. We also studied the invariant mass distributions for the two-body subsystems to 
search for new resonances. The analysis is performed on a data sample containing (657 ± 
9) X 10^ BB pairs, collected with the Belle detector at the KEKB asymmetric-energy e^e~ 
collider [4] that operates at the T(4S') resonance. The production of B^B~ and B^B pairs 
is assumed to be equal. 



b 



B+ 




FIG. 1: Diagrams for the decays (a) B+ Dy' K+t:+ and (b) B+ Dr'^D^.D^ K+tt-. 



The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) composed of CsI(Tl) crystals, located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron fiux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [5]. Two inner detector configurations 
were used. A 2.0 cm beam pipe and a 3-layer silicon vertex detector were used for the first 
sample of 152 x lO^BB pairs, while a 1.5 cm beam pipe, a 4-layer silicon detector and a 
small-cell inner drift chamber were used to record the remaining 505 x 10^ BB pairs [6]. 

Charged tracks are required to have a distance of the closest approach to the interaction 
point less than 5 cm in the beam direction (along the z-axis) and less than 5 mm in the 
transverse (r — 0) plane. In addition, we only select charged tracks that have transverse 
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FIG. 2: Distributions of AE, Mbc and M(L>^) for (a) B+ D-{-^ (I)-k-)K+tt+ , (b) B+ 
D-{^ K*^K-)K+7r+ and (c) B+ D-{-^ K'^K-)K+tt+ decays. Each distribution, AE, Mbc 
or M{Ds) includes a selection on the signal region of the remaining two. The red dashed curves 
show the results of the overall fit described in the text, the blue dashed curves correspond to the 
signal components and the green dashed curves indicate the fitted background for M{Ds)- 



momenta larger than 100 MeV/c. 

To identify charged hadrons, we combine information from the CDC, ACC and TOF into 
pion, kaon and proton likelihood variables ^^k and Cp. For kaon candidates we then 
require the likelihood ratio Ck/tt = r '~lr to be larger than 0.6. We also apply the proton 
veto condition: Lpji^ < 0.95. Pions are selected from tracks with low kaon probabilities 
satisfying a likelihood ratio condition Cx/n < 0.6 together with a proton veto Cp/K < 0.95. 
In addition, we reject all charged tracks consistent with the electron or muon hypothesis. 
The above selection results in a typical kaon (pion) identification efficiency ranging from 
92% to 97% (94% to 98%) for various decay modes, while 2% to 15% of kaon candidates are 
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misidentified pions and 4% to 8% of pion candidates arc misidcntificd kaons. 

The candidates are reconstructed in three final states: K^)!^^ ^ iir*(892)°(— >■ 

K-t:+)K+ and K^s{-^ 71+7t-)K+. We accept K+R- {K-71+) pairs as (p_(K^{892f) candi- 
dates if their invariant mass is within 10 (100) MeV/c^ of the nominal (A'* (892)") mass [7]. 
This requirement corresponds to ±2.5(7 in all cases. Candidate Kg mesons are selected by 
combining oppositely charged particles with an invariant mass not differing by more than 
6 MeV/c^ from the nominal Kg mass. In addition, the vertex of these tt+tt" pairs must 
be displaced from the interaction point by at least 5 mm. Photons used for D* — > Dg^ 
reconstruction are accepted if their energies exceed 100 MeV in the laboratory frame. No 
selection requirements arc imposed on the ui*"^ mass at this stage. 

A B meson is reconstructed by combining the Dg candidate with an identified kaon and 
pion and by applying a loose requirement on the quality {xb) of the vertex fit to the K, n and 

Dg trajectories, where the Dg mass is constrained to the world average value [7]. The signal B 

(*) 

meson decays are identified by three kinematic variables, the Dg invariant mas s, the energy 
difference, AE = Eb — Ei^csun- and the beam-energy-constrained mass, Afbc = \/ -^bcam ~ Pb- 
Here Eb and ps arc the reconstructed energy and momentum of the B candidate, and i?beam 
is the run-dependent beam energy, all are calculated in the center-of-mass (CM) frame. For 
further analysis we retain events in the candidate region defined as: 1.91 GeV/c^ < M{Dg) < 
2.03 GeV/c^ (2.06 GeV/c^ < M{D*) < 2.16 GeV/c^), 5.2 GeV/c^ < Mbc < 5.3 GeV/c^ 
and —0.08 GeV < AE < 0.2 GcV. The lower bound in AE for candidate events is chosen 
to exclude a possible background from B DgX decays with higher nniltiplicitics. From 
GEANT [8] based Monte Carlo (MC) simulation, we deduce that the signal peaks in a 
signal box are defined by the requirements: 1.9532 GeV/c^ < M{Dg) < 1.9832 GeV/c^ 
(2.092 GeV/c2 < M{Dl) < 2.132 GeV/c^), 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ and \AE\ < 
0.03 GeV. Based on MC simulation, the region 2.88 GeV/c^ < M{K+K-7r+7r-) < 3.18 
GeV/c^ is excluded to remove background from B^ — > ((cc) — >■ K~^K~7i~^7r~)K~^ decays, 
where (cc) are charmonium states such as the J/ip or rjc- For 5+ — > DfD^{^ X+tt") decays 
the events in the candidate region are required to have the K'^tt^ invariant mass within a 
15 MeV/c^ (3(t) interval of the nominal mass. 

We find that for B+ D-K+n+ {B+ D*-K+n+) decays at most 11% (29%) of 
events have more than one B candidate. In such cases we select the B candidate with the 
smallest value of x%- Moreover, when there are at least two combinations with the same x'b 
value, the one containing a kaon - originating directly from the B decay - with the highest 
likelihood ratio Ck/w is selected. For B D*K7r decays, wc further choose the combination 
that minimizes the quantity \M{D*) - M{Dg) - 143.8 MeV/c^]. _ 

We exploit the event topology to discriminate between spherical BB events and the 
dominant background from jet-like continuum events, e'*"e~ Qq {q — u, d, s, c). We 
use the event shape variable R2 defined as the ratio of the second and zeroth Fox- Wolfram 
moments [9] and require that R2 be less than 0.4. 

The signal yields are extracted using unbinned extended maximum-likelihood fits to the 
(AE , M-bc, M {ni^'' )) distributions of the selected candidate events. The hkehhood function 
is given by: 

C = ^(Ns ^ TV.J-'e— - n {j^n + . (1) 

where i is the event identifier, N is the total number of events in the fit and Ns{Nb) is 
the number of signal and background events, respectively. We use Gaussian functions to 
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FIG. 3: Distributions of AE, Mbc and M{D*) for (a) ^ I?r(^ </>7r-)i^+7r+, (b) 
-D*-(^ K*°K-)K+7r+ and (c) -B+ ^ I)r(^ A'g/i:-)i^+7r+ decays. For each distribution, AE, 
-/Vfbc and M{D*) we select on the signal region of the remaining two. The red dashed curves show 
the results of the overall fit described in the text, the blue dashed curves correspond to signal 
components and the green dashed curves indicate the fitted background for M(D*). 



parameterize the signal probability density function in AE and Mbc and a double Gaussian 
function with a common mean for the M(.Di*^) distribution: 



X /^(.) GiM^D^y, + GiM^D^y, 

where AE,mB,m {,),aAE,o'Mi,^, f^(,),cr^^},) and 17^%) are fit parameters. The latter three, 
which describe the signal shape corresponding to the mDi*^ distributions are fixed to the 
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values obtained from the fit to the S+ — > D^*^'^^ control channels. In addition, we use the 
B+ D*+^ data samples to fix the signal widths for AE and Mbc for the B+ D*-K+tt+ 
decays. 

The background is parameterized with a second-order polynomial (p2) in the Ai5 distribu- 
tion. For the Afbc background distribution we choose a parameterization t hat was first used 
by the ARGUS collaboration [10], f{M^eX) oc M^,^l - {M^J Ey,,^Ye-<^^-^''^^'''^^--'^"\ 

where is a fit parameter. Finally, the M(Ds*'') background distribution is described by the 
sum of a double Gaussian function and a second-order polynomial: 

n =P2(AE^ wo,w,,W2) X /(M^,; C)x (3) 



X 



The values of the variables Wq, Wi, t(72, C, "^Oi "^ii "^2 are determined in the fit, whereas the 



arc fixed to the values resulting from the fits to the appropriate control channels. Figures 2 
and 3 show the distributions of AE, Mbc and M{D't'') together with the fits described above. 

For decays containing a K*{d>92)^ meson a small correction was apphed to the signal yields 
obtained from the fit. The /s:*(892)° mass sidebands (0.746-0.796) GeV/c^ and (0.996-1.046) 
GeV/c^ were fitted and a significant background contributing to the signal yields was found 
for the B+ D-{^ K*^K-)K+7r+ and B+ D+{^ K*^K-)DO channels. Final signal 
yields were obtained by subtracting these contributions from the nominal fit values. 

The signal yields tog ether with stati stical significances are fisted in Table I. The sig- 
nificance is defined as a/— 21n(£o / -Cmax) , where £inax (-^o) denotes the maximum likelihood 
with the signal yield at its nominal value (fixed to zero). 

The reconstruction efficiencies, determined using MC samples of e"'"e~ — > T(45') — > B^B~ 
decays, are listed in Table I. This table also contains the values obtained for the branching 
fi-actions of the decays 5+ Di*^~K+ 7r+ and B+ Di*^^W. The last error (Table I) 
is due to uncertainties in the branching fractions for the decays of intermediate particles, 
predominantly the oi*"^ [7]. The systematic uncertainties are evaluated only for the three- 
body B~^ — > Di*^~ K~^7r~^ decays. We find branching fractions for the B~^ — > Di*^~^D^ control 
samples in agreement with world averages [7] . 

Systematic uncertainties are fisted in Table II. The contribution (f ) due to the selection 
procedure is dominated by the R2 requirement. This uncertainty is estimated conservatively 
as the maximum variation of the efficiency-corrected signal yield, when the R2 selection value 
is varied over a wide range (values between 0.25 and 0.55). The uncertainty (g) due to the 
fit range is determined by varying the candidate region. To evaluate the contribution (h) 
we repeat the fits varying the shape parameters by ilex. The uncertainty (i) is estimated as 
the statistical error in the selection efficiency, increased conservatively by a factor obtained 
from the difference between the value of the branching fraction for the appropriate control 
channel and the generated branching fraction. The overall systematic error is obtained by 
summing these contributions in quadrature. 

The average branching fractions for the decays B^ — > D~K^7r^ and B~^ D*~K^7r~^ 
are determined from a simultaneous fit to the data containing events from all three Dg 
decay modes. Here, the systematic uncertainties are calculated as in the individual channels 
(Table II). 
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TABLE I: Signal yields, reconstruction efficiencies, branching fractions and statistical significances 
for B+ Di*^'K+TT+ and B+ L>i*^+K+7r- decays. 



Decay 



Signal Efficiency Statistical 



Branching 



[%] Sigiiif. [a] fracliou [(IQ-^)] 



B+ - 






(f)Tr )K~^Tr~^ 




11.25 ±0.11 


31.5 


i.9o; 


f 0.12 +0.18 1 n OQ 
-0.12 -0.18 ^ ^-^^ 


B+ - 


-D- 


(- 


K*^K-)K+TT+ 




8.55 ±0.10 


26.5 


1.93; 


f0.17 +0.19 1 n on 
-0.16 -0.19 ^ "-"^^ 


B+ - 


-D- 


(- 






12.82 ±0.19 


20.4 


2.06; 


f0.19 +0.25 _i_n TO 
-0.18 -0.26 ^ ^■^'^ 


B+ - 


-Dt 


-(- 






3.00 ± 0.06 


11.0 


1.46; 


f 0.23 +0.18 1 A 99 
-0.21 -0.19 ^ ^-^^ 


B+ - 








O-*^- ' -9.8 


2.65 ± 0.06 


9.3 


1.45; 


f 0.25 +0.18 1 n 99 
-0.23 -0.19 ^ ^-^^ 


B+ - 




-(- 




or; 7+7.7 


3.67 ±0.11 


8.0 


1.53; 


f 0.33 +0.24 1 n no 
-0.30 -0.22 ^-^^ 



B+ Df{-^ 07r+)L>o 
B+ D+{-^ 'K*^K+)W 
B+ ^ D+{^ K^K+)W 

B+ D*+{-^ (l)TT+)W 

B+ ~ 
B+ 



s 

D*+ 



597.4+25 13.03 ± 0.12 

512.61^1^ 9.21 ±0.10 

294.5li^:^ 14.22 ± 0.20 

150.2tlli 3.97 ±0.07 

151.911^;.^ 3.09 ±0.06 

^^■^-iii 4.40 ±0.12 



56.8 
53.3 
38.9 
19.0 
20.8 
15.0 



82.3ll|3^ ± 12.50 
83.80lt?! ± 12.94 
78.61^4:16 ± 4.86 
72.151^:^2 ± 10.98 

78.681^:^^ ± 12.16 
87.27liJ:^^ ± 5.43 



TABLE II: Systematic uncertainties on the branching fractions for 5+ —>■ Dg K'^tt'^ decay 
modes, given in percent. 



Source 


D 


g final state 


D*~ 

^ s 


final state 




(pTr~ 


K*^K- 




^7r^7 K 






(a) Tracking 


5 


5 


5 


5 


5 


5 


(b) Hadron identification 


5 


5 


5 


5 


5 


5 


(c) Kg reconstruction 






4.5 






4.5 


(d) Photon reconstruction 








5 


5 


5 


(e) Uncertainty in N{BB) 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


(f) Selection procedure 


3.7 


3.7 


3.7 


3.7 


3.7 


3.7 


(g) Size of candidate region 


0.6 


1.1 


1.2 


0.3 


0.2 


1.5 


(h) Signal shape 


+1.2 

-1.3 


+3.4 

-3.4 


+5.6 
-6.3 


+5.6 
-6.7 


+6.9 
-7.4 


+9.7 
-7.8 


(i) MC statistics 


5.1 


4.1 


5.2 


5.4 


4.5 


5.8 


Total 


+9.6 
-9.7 


+9.8 
-9.8 


+12.1 
-12.4 


+12.3 
-12.8 


+12.6 
-12.9 


+15.5 
-14.4 



In summary, the following branching fractions are determined: 

B{B+ ^ D:K+7t^) = (1.94l°:°^(stot) t'o^isyst) ± 0.17(23,,*)) x 10"^ (4) 

B{B+ ^ L>rx+7r+) = (1.47l°:l^(stoi) t^^^^{syst) ± 0.13(Bi„t)) x 10"^ (5) 

These branching fractions are compatible with the values reported by the BaBar collabora- 
tion [2]. 
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The invariant mass distributions of the Di*''~ K'^ subsystem are incompatible with those 
expected for three-body phase space production and exhibit strong enhancements around 
2.7 GeV/c^ (see Fig. 4). These features may be explained by the production of charm 
resonances with masses below Di*^ K+ threshold [3]. 
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FIG. 4: The invariant mass distributions of (a) D^K^ for the decay D^K^ir^ and (b) 

D*^ K^TT^ corresponding to the signal regions described in the text. The 
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histograms show the background contributions corresponding to AE > 0.06 GeV. 
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